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Abstract 

Issues of morphology, nucleation and growth of Ge cluster arrays deposited by ultrahigh vacuum molecular beam 
epitaxy on the Si(OOl) surface are considered. Difference in nucleation of quantum dots during Ge deposition at 
low (<600°C) and high (>600°C) temperatures is studied by high resolution scanning tunneling microscopy. 
The atomic models of growth of both species of Ge huts — pyramids and wedges — are proposed. The growth 
cycle of Ge QD arrays at low temperatures is explored. A problem of lowering of the array formation 
temperature is discussed with the focus on CMOS compatibility of the entire process; a special attention is paid 
upon approaches to reduction of treatment temperature during the Si(OOl) surface pre-growth cleaning, which is 
at once a key and the highest-temperature phase of the Ge/Si(001) quantum dot dense array formation process. 
The temperature of the Si clean surface preparation, the final high-temperature step of which is, as a rule, 
carried out directly in the MBE chamber just before the structure deposition, determines the compatibility of 
formation process of Ge-QD-array based devices with the CMOS manufacturing cycle. Silicon surface 
hydrogenation at the final stage of its wet chemical etching during the preliminary cleaning is proposed as a 
possible way of efficient reduction of the Si wafer pre-growth annealing temperature. 



Introduction: Background and problem statement 



Heteroepitaxial structures Ge/Si and SiGe/Si are among the most promising materials of modern 
nanoelectronics and nanophotonics pHT2]. Lately industry has developed numerous radiofrequency devices 
on the basis of SiGe/Si structures with bands wider than 100 GHz, which already compete with 
Ga As-based components. Except that, SiGe-based technology has allowed one to approach to development 
of the most important elements of single-crystalline integrated microphotonics — laser diodes and detectors 
for fiberoptic communications; SiGe-based waveguides are already available. So, the forecasted forthcoming 
breakthrough would enable the solution of two important problems — (i) development of monolithic VLSI 
circuits for fiberoptic telecommunications and (ii) replacement of electronic data buses by optical ones. 
Additionally, encouraging results have been recently obtained in development of emitting THz and 
mid-infrared devices based on SiGe/Si- heterostructures. And finally, application of Ge/Si and SiGe/Si 
hetersructures might enable a breaking progress in IR imaging technology opening a way to creation of 
multispectral photodetector arrays integrated with readout circuitry on a single-crystalline chip. 

Dense arrays of Ge quantum dots (QD) are of importance to all practically significant applications in 
optoelectronics and microelectronics. QD array is usually referred to as dense array if interactions among 
adjacent clusters play an important role [13], i.e., a tunnel coupling between Ge clusters arise [H]; such 
arrays should be considered as a whole in terms of behaviour of current carrier and transport 
properties [10]. Ge/Si hetrostructure can include both isolated QD arrays, i.e., arrays which do not 
mutually interact and separated by a thick enough (> 30 nm) Si layer, and superlattices of QD arrays, in 
which arrays are separated by thin Si layers and which, like SiGe/Si superlattices with quantum wells, 
represent a single coupled system. In addition, it is known that, if distances between layers of Ge clusters 
are as small as unities of nanometers in the direction of the structure growth, ordering of Ge clusters is 
observed in this direction [15 ; they form chains which can be composed by dozens of clusters if a number 
of layers is large [16 . Like atoms, Ge clusters form a sort of molecules in which electron density 
redistributes among clusters depending on distances between them (thicknesses of Si buffer layers), as if 
changing a type of chemical bond from covalent to ionic. This phenomenon opens a wide perspective to 
designing heterostructures with various optical and electrical properties. 

Recently, an interest of researchers has been attracter by heterostructures with ordered arrays of quantum 
dots, in particular, by ordered arrays of Ge nanoclusters in the Si matrix [T7]- An idea of formation of a 



2 



QD array, which would combine advantages of a single QD with benefits of a dense array, seems to be the 
most promising.^ Controllable ordering of clusters in all three directions would enable creation of a 
volumetric crystal in which QDs play a role of atoms. In such artificial crystal, unique opportunities 
appear to designing wave functions of carriers by filling corresponding quantum states of QDs by electrons 
or holes (like it is the case for 5, p and d atomic configurations). As opposed to impurity states in 
semiconductors, dense 3D array of QDs would be an ensemble of multicharged centres in which an essential 
role would be played by the Coulomb potential. A concept of the QD crystal, which is considered as a 3D 
lattice of artificial atoms, implies a new material with spatial ordering on the scale comparable with the 
de Broglie wave-length for electrons. Non- locality of the quantum- mechanical bonding together with the 
Coulomb interaction of carriers localized in close QDs may result in new optical and electronic properties 
arising from the collective nature of electronic states. In contrast to stochastically located QDs, in this case 
these properties would not be averaged over components of a crystal. Main properties of such ensemble 
would reproduce peculiarities inherent to ordinary solids, such as appearance of two-dimentional or 
three- dimentional minibands, separated by minigaps, in lieu of localized quantum states intrinsic to 
separate QDs. It is necessary for QD crystal that QDs would be ordered to precise periodicity, the sizes of 
QDs would equal, and distances between QDs would be small enough for wave functions to overlap. 

Such QD crystals would be very prospective for application in nanoelectronics, spintronics and, likely, in 
quantum computing, as well as in devices of silicon optoelectronics such as highly efficient sources and 
detectors of infrared and terahertz emission enabling integration to silicon VLSI circuitry. 

Main restriction for use of such Ge/Si heterostructures with dence arrays of self arranged Ge QDs are 
associated with the spread sizes of Ge clusters and their tend to disordering on the surface. Both these 
factors cause tailing of a discrete spectrum. Additional difficulty is the necessity for all technological steps 
to be embedded into VLSI manufacturing process or, in other words, meet requirements of CMOS 
compatibility. 

To be able to solve the above ambitious task, a deep knowledge of physical processes on silicon surface 
during its preparation and in germanium and silicon films during the heterostructure formation is strongly 
required. This article represents some results of our resent investigations in this direction. 
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Methods 
Equipment 

The experiments were carried out using an integrated ultrahigh vacuum instrument p!8H2Q] built on the 
basis of the Riber SSC 2 surface science center with the EVA 32 molecular-beam epitaxy (MBE) chamber 
equipped with the Staib Instruments RH20 diffractometer of reflected high-energy electrons (RHEED) and 
connected through a transfer line to the STM GPI-300 ultrahigh vacuum scanning tunnelling 
microscope [2TH23].^ A preliminary annealing and outgassing chamber is also available in the instrument. 

The pressure of about 5 x 10~^ Torr was kept in the preliminary annealing chamber. The MBE chamber 
was evacuated down to about 10~^^ Torr before processes; the pressure grew to nearly 2 x 10~^ Torr at 
most during the sample surface deoxidization process and 10~^ Torr during Ge or Si deposition. The 
residual gas pressure did not exceed 10~^^ Torr in the STM chamber. 

The instrument enables the STM study of samples at any stage of Si surface preparation and MBE growth. 
The samples can be consecutively moved into the STM chamber for the analysis and back into the MBE 
vessel for further treatments or Ge, Si or SiGe deposition as many times as required never leaving the UHV 
ambient and preserving the required cleanness for MBE growth and STM investigations with atomic 
resolution. RHEED experiments can be carried out in situ^ i.e., directly in the MBE chamber during a 
process [T9] . 

Sources with the electron beam evaporation were used for Ge or Si deposition. The deposition rate and 
coverage were measured using the Inficon Leybold-Heraeus XTC 751-001-Gl film thickness monitor 
equipped with the graduated in-advance quartz sensors installed in the MBE chamber. Tantalum radiators 
were used for sample heating from the rear side in both preliminary annealing and MBE chambers. The 
temperature was monitored with chromel-alumel and tungsten-rhenium thermocouples of the heaters in 
the preliminary annealing and MBE chambers, respectively. The thermocouples were mounted in vacuum 
near the rear side of the samples and in situ graduated beforehand against the IMPAC IS 12- Si pyrometer 
which measured the sample temperature through chamber windows. The temperature distribution 
uniformity over a surface was also investigated in advance; the deviations from mean values were found to 
be within ±3°C for the half-radius areas around the centers of 2'' wafers over the whole temperature 
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interval applied in this study. 

The composition of residual atmosphere in the MBE camber was monitored using the SRS RGA-200 
residual gas analyzer before and during the process. 

The STM tip was ex situ made of the tungsten wire and cleaned by ion bombardment [2T in a special 
UHV chamber connected to the STM one. 

In this work, the images were obtained in the constant tunnelling current (/t) mode at the room 
temperature. The STM tip was zero-biased while the sample was positively or negatively biased (t/g) when 
scanned in empty- or filled-states imaging mode. 

Original firmware [2TH23] was used for data acquisition; the STM images were processed afterward using 
the WSxM software [25]. 

Sample preparation procedures 

Preparation of samples with deposited Ge layers 

Initial samples for STM were 8x8 mm^ squares cut from the specially treated commercial boron-doped 
Czochralski-grown (CZ) Si(lOO) wafers (p-type, p = 12 l^cm). After washing and chemical treatment 
following the standard procedure described elsewhere [26 (which included washing in ethanol, etching in 
the mixture of HNO3 and HF and rinsing in the deionized water [27 ), the silicon substrates were mounted 
on the molybdenum STM holders and inflexibly clamped with the tantalum fasteners. The STM holders 
were placed in the holders for MBE made of molybdenum with tantalum inserts. Then, the substrates were 
loaded into the airlock and transferred into the preliminary annealing chamber where they were outgassed 
at the temperature of around 565°C for more than 6h. After that, the substrates were moved for final 
treatment and Ge deposition into the MBE chamber where they were subjected to two-stages anneahng 
during heating with stoppages at 600°C for 5 min and at 800°C for 3 min [18] . The final annealing at the 
temperature greater than 900°C was carried out for nearly 2.5 min with the maximum temperature of 
about 925°C (1.5 min). Then, the temperature was rapidly lowered to about 750°C. The rate of the further 
cooling was around 0.4°C/s that corresponded to the 'quenching' mode applied in [19 . The surfaces of the 
silicon substrates were completely purified of the oxide film as a result of this treatment [T9l[28l[29] . 

Ge was deposited directly on the deoxidized Si (001) surface. The deposition rate was varied from about 0.1 
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to 0.15 A/s; the effective Ge film thickness {hce) was varied from 3 to 18 A for different samples. The 
substrate temperature during Ge deposition (Tg^) was 360°C or 530°C for the low-temperature mode and 
600 or 650°C for the high-temperature mode. The rate of the sample cooling down to the room 
temperature was approximately 0.4°C/s after the deposition. 

After cooling, the prepared samples with Ge layers were moved for analysis into the STM chamber. 
Preparation of samples for Si(OOl) surface analysis 

Wafers for Si (001) surface analysis by STM and RHEED were the same as for Ge MBE. Initially, the 
specimens were chemically etched in the RCA etchant [30 and processed to form a surface terminated by 
hydrogen atoms (Si:H). The hydrogenated Si:H samples were prepared by etching in solutions containing 
HF at the final stage of the RCA process [31 . We used the following solutions with pH = 2, 4 or 7: a dilute 
HF solution (5% or 0.5%), buffered NH4F + HF or NH4F solutions. After that, the Si:H samples were 
pre-treated for 2 hours at the temperature of ~300°C and the pressure of less than 5 x 10~^^ Torr in the 
MBE chamber. The second phase of the thermal treatment was conducted at the temperatures of 800, 650, 
610, 570, 550, 530 or 470°C. Duration of this phase was chosen to a form of the RHEED pattern. The 
samples were quenched after heat treatments at the rate of ~0.4 °C/s [19] . 

Results and Discussion 
Nucleation 

Hut nucleation at low and high temperatures 

Nucleation of Ge clusters at low temperatures has been an issue of numerous experimental and theoretical 
investigations for a number of years (see a brief review section in article [H]). Recently, we have described 
two characteristic formations composed by epitaxially oriented Ge dimer pairs and chains of four dimers on 
the wetting layer (WL) patches which were interpreted by us as two types of hut nuclei: an individual type 
for each species of huts — pyramids or wedges (Figure [Bi-c) [T8l[20l[32] . These nuclei are always observed to 
arise on sufficiently large WL patches: there must be enough room for a nucleus on a single patch; a 
nucleus cannot be housed on more than one patch [32 . 

Both types of the hut nuclei appeared to arise at the same WL M x N patch thickness [20l[33], hence, at 
the same WL stress to relieve it. Therefore, they appear at the same strain energy (and with equal 
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likelihoods, see Refs. [T8l[32]). This means that they are degenerate by the formation energy: if they had 
different formation energies they would appear at different WL thicknesses; the first of the types of huts, 
which nucleates on the surface, releases the stress; the second one never appears therefore. Hence, they can 
occur only simultaneously and their formation energies can be only equal. 

Presently we have no satisfactory explanation of this phenomenon and can only propose a very preliminary 
interpretation of the observed simultaneous appearance of the two kinds of nuclei on WL. The explanation 
is based on modeling of Ge cluster formation energy performed in Ref. [34 . Brehm et al. [S^ have explored 
Ge island nucleation during MBE at much higher temperatures than those applied in this work, therefore 
theoretical results of Ref. [34 describe the experimental data obtained for the case of the high-temperature 
growth mode, which differs considerably from the low-temperature one [18 . However, the modeling could 
also apply for the low-temperature growth. According to Ref. [34 , flat Ge islands — in our case, nuclei and 
small huts — likely occur on WL because of an energy benefit which arises in exposing the compressed 
{105} facets, rather than in relaxing the volumetric elastic energy, as it takes place in the usual 
Stranski-Krastanov mechanism. At low temperatures, this effect may stabilize clusters, however preventing 
their further ripening (this agrees with our observations presented recently in Ref. [18 ). If this is the case, 
the actual volumetric and structural form of clusters likely do not impact very much in their formation 
energy.^ 

As distinct from the low-temperature mode, Ge cluster nucleation at high temperatures may go on in two 
ways. The first way is similar to the process of hut nucleation at low temperatures. Pyramids were 
observed to nucleate in such a way. Figure [l]i,e illustrates this statement: the pyramid nuclei, absolutely 
the same as those observed in the samples grown at low temperature, are seen on the WL patches in the 
images of the samples obtained at Tgr = 650°C. Their density was small, and they were mainly situated in 
the vicinity to large mature pyramids, which arise at early stages of Ge deposition and have much greater 
sizes than huts formed at low temperatures at the same values of Hqq [35j. The WL surface mainly 
consisted of monoatomic steps and narrow terraces in these ares (Figure [T]i). 

The second way, somewhat resembling the process described by Goldfarb et al. [36 for the case of the 
gas-source MBE (and thick hydrogenated WL), is illustrated by Figure O At small values of /iQe, regions 
containing excess of Ge atoms were observed on the surface. Usually, they were not resolved as structured 
formations and resembled shapeless heaps of Ge (Figure [2^). Pits usually accompanied them. Heap density 



7 



was about 10^ cm ^. Some of heaps had started to form the {105} facets during Ge deposition (Figure [2]3). 

Stoppage of Ge deposition and subsequent anneahng at Tgr resulted in formation of volumetric structures 
partially or even completely faceted by {105} planes, transforming 'heaps' to some similarities of huts 
(Figure [2fc,d,e). 

We have never observed such process at low temperatures of growth and suppose it to be inherent only to 
the high-temperature array formation mode. 

Array nucleation and growth outset 

Since the pioneering work by Mo et al. |37], it has been known that deposition of Ge on Si(OOl) beyond 
3 ML (1 ML ^1.4 A) leads to formation of huts [37H39] on WL with high number density (> 10^° cm-^ 
Refs. [T8l[20lll0]). Some later the value of Ge coverage, at which 3D clusters emerged, was confirmed by 
Iwawaki et al. [JT] who, in the course of a comprehensive STM study of the low-temperature epitaxial 
growth of Ge on Si(OOl) j4TH44] . directly observed appearance of minute (a few ML height) 3D Ge islands 
at 300°C on (M x A/')-patched WL; deposition of 4 ML of Ge resulted in formation of a dense array of small 
huts. Various values of Ge coverage, at which the transition from 2D to 3D growth occurs, are presented in 
the literature. For example, an abrupt increase in hut density at the coverage of 3.16 ML was detected for 
Ge deposition at 300°C and 0.06ML/min [40 . A detailed phase diagram of the Ge film on Si(OOl) derived 
from experiments carried out by recording RHEED gave the coverages corresponding to the "2D-to-hut" 
transition from ~2.5 to ~3ML for the growth temperature interval from 300 to 400°C (and different values 
for different temperatures) [45]. Photoluminescence study of Ge huts deposited at the temperature of 
360°C showed that evolution from "quantum-well-like" (attributed to WL) to "quantum-dot-like" 
(attributed to Ge huts) emission occurred at a coverage of '^4.7 ML in PL spectra obtained at 8K [46] . 
Hut formation studied by high resolution low-energy electron diffraction and surface-stress-induced optical 
deflection evidenced that at deposition temperature of 500°C hut formation suddenly set in at a coverage of 
3.5 ML [47 . And finally, for theoretical studies the WL thickness and consequently the hut formation 
coverage is usually assumed to equal 3 ML [48 . As it is seen from the above examples, there is no 
unambiguous information presently about the coverage at which huts arise or, more accurately, about the 
thickness of the WL M x N patch on which a cluster nucleate during Ge deposition. STM studies show the 
WL thickness to equal 3 ML only on the average: M x N patches have slightly different thicknesses 
(dz 1 ML) around this value [T8l[20l[32j|lTll39] . In this section, we determine by means of high resolution 
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STM an accurate value of the M x N patch height at which hut nucleate at 360°C. 

Figure [3^1 demonstrates a typical STM microg raph of the (M x Ar)-patched WL {hce = 4.4 A, -3.1 ML). 
This image does not demonstrate any feature which might be recognized as a hut nucleus (Figure [T^-c) [32] . 
Such features first arise at the coverages of —5 A: they are clearly seen in Figures [Sb-d, which demonstrate 
a moment when the array have just nucleated {hce = 5.1 A, —3.6 ML). However, we succeeded to find 
minute pyramid and wedge at this hce (Figure [3]i) — both as small as 2 ML over the patch surface (we 
measure cluster heighs from patch tops) — which indicate that hut nucleation had started a little earlier. 

It can be concluded from these observations that hut arrays nucleate at a coverage of -5.1 A (-3.6 ML) 
when approximately a half of patches are as thick as 4 ML. We can suppose then that huts nucleate on 
those patches whose thickness have reached (or even have exceeded) 4 ML. 

Morphology 

Wetting layer reconstruction 

Evolution of WL patches during MBE is illustrated by Fig. HJ In full agreement with the data of Ref. [?T], 
both c(4 X 2) and p{2 x 2) reconstructions are observed on tops of the M x N patches in all images except 
for the image given in Figured^ {hce = 4.4 A) in which only the c(4 x 2) structure is recognized. A 
magnified image of the p{2 x 2) structure illustrating its characteristic zig-zagged shape and resolving 
separate upper atoms of buckled Ge dimers is given in Figure |3]i. 

Formation of a hut nucleus on a patch reconstructs its surface; a new formation changes the structure of 
the topmost layer to that specific for a particular type of nuclei, in the present case, to the structure of the 
pyramidal hut nucleus (Figure [4^). However the residual c(4 x 2) structure still remains on the lower 
terrace of the patch. At the same time, the p{2 x 2) structure stays on the top of the adjacent patch [33j. 

We can conclude now that c(4 x 2) and p{2 x 2) surface structures occurring on the M x N patches should 
be energetically degenerate. In addition, the above observation rises a question if there is some connection 
between the form of a patch top reconstruction and a species of hut which could nucleate on it or, in other 
words, whether the patch top reconstruction controls hut nucleation and determines its species. 
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Growth and structure of pyramids and wedges 

Let us consider possible scenarios of hut growth after nucleation on WL. Earher [18l|20l|49] , we have 
aheady proposed structural models of both species of huts and briefly discussed processes and atomic 
models of their formation giving a few drawings with identical apexes as examples and allowing the readers 
to construct the missing structural schemes. However, crystallography allows one to arrive to two different 
solutions for wedge-like huts, and additional empirical knowledge and STM data is required to discriminate 
between them. Both solutions are given in Figure [5l The first scenario of growth assumes uniform addition 
of Ge atoms to all four facets of huts (follow a series number I in Figure [5j). In this case wedge-like huts 
have different ridge structures (the ridge width and location of atoms on it) depending on cluster height. 
The initial ridge structure, which form on top of 2-ML wedge reconstructing the nucleus [20l[32], should 
then occur on the ridge every 5 ML over the nucleus, i.e., only the wedges of 2, 6, 11, etc. ML height over 
WL can have the same ridge structure. This contradicts our observations according to which the structure 
of hut apexes always remains and depends on only the hut species. Therefore, we are made to come to a 
different solution in which Ge atoms are added to facets non-uniformly (see series II for wedges and 
pyramids in Figure [5]). In this scenario, the structures of both apexes of huts are independent of cluster 
heights, that agrees with experimental observations. 

Complete cycle of Ge QD array growth at low temperatures 

The STM images of the surfaces of the germanium layers grown at Tgr = 360°C with various Hqq values are 
shown in Figure [6] where the evolution of the Ge layer on the Si(OOl) surface in the process of 
low-temperature MBE is seen. Hut clusters on the Ge surface have not nucleated at Hqq — 4.4 A, and the 
STM image in Figure [6^ exhibits only the well-known structure of the wetting layer with the c(4 x 2) or 
p{2 X 2) reconstruction inside M x N blocks [20l[32]. The array nucleates at hce ~ 5 A (Figure [3|) but 3D 
huts mainly form at higher coverages [20[[33]. Hut arrays initially evolve with increasing hce by concurrent 
growth of available clusters and nucleation of new ones resulting in progressive rise of hut number density. 
Huts are clearly seen in Figure [GJd for hce = 6 A; their density and sizes increase; the number density of 
huts reaches maximum at hce = 8 A (Figure [Bt); clusters with various sizes — completely formed clusters, 
recently nucleated small clusters, and nuclei with a height of 1 ML over the Ge WL — are simultaneously 
present on the surface j20l|32l[33] . This array is very inhomogeneous both in the sizes of the clusters and in 
composition; it includes regular pyramidal and elongated wedge-shaped clusters, but wedge-shaped clusters 
with a large spread in the lengths dominate [18]. The array is most homogeneous at hce = 10 A 
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(Figure [Gji) [50], clusters cover almost the entire surface of the wetting layer, the fraction of small clusters 
decreases noticeably, and large clusters begin to coalesce. At hce = 14 A, most clusters coalesce near their 
bases (Figure [6^,f), and the free wetting layer almost disappears from the field of view of STM, but the 
array consists of individual clusters. At hce = 15 A, the coalescence of clusters continues and a transition 
to the growth of a two-dimensional film of nanocrystalline germanium begins (Figure [6^). Nevertheless, the 
hut nucleation continues on small lawns of WL rarely preserved, surrounded by large huts, even at as high 
coverages as 15 A, when virtually total coalescence of the mature huts have already happened [20j. Finally, 
at hce = 18 A, it is seen that the array of Ge clusters disappears and although the roughness of the surface 
is still pronounced, the Ge layer grows as a continuous nanocrystalline film (Figure [6]i). A chaotic 
conglomeration of faceted hillocks and pits composes the film; steep facets appear around the pits 
(Figure [6j). However, Ge WL (M x A^)-patched structure is clearly resolved on the bottom of pits and WL 
lawns (Figure [6|,k). WL appears to be a very stable formation. 

The density of the wedge-shaped clusters increases when hce increases up to 8 A and, then, decreases 
slowly, whereas the density of the pyramid-shaped clusters decreases exponentially in the process of growth 
of the array jl8[l2Uj. The total density of the clusters is about 3.5 x 10^\ 5.8 x 10^\ 5.1 x 10^\ and 
2.3 X 10^^ cm~^ at hce = 6, 8, 10, and 14 A, respectively. From the capacitance- voltage characteristics of 
the samples with Ge/Si(001) heterostructures, the surface densities of holes in them were earlier estimated 
as 3.4 X 10^^, 7.0 x 10^^, and 1.7 x 10^^ cm~^ for hce = 6, 10, and 14 A, respectively. These values almost 
coincide with the densities of the Ge clusters in arrays [5T1[52]. Notice also, that very high terahertz 
conductivity was observed by Zhukova et al. [52 in the samples with hce = 8, 9, 10, and 14 A which 
drastically decreased for hce = 18 A and was not detected at all at 6 A and lower values of hce- 

CMOS compatibility 

CMOS compatibility of technological processes based on Ge/Si heteroepitaxy imposes a hard constraint on 
conditions of all the phases of the heterostructure formation including Si wafer thermal cleaning and 
surface preparation to epitaxial growth. Formation of a device structure with QD arrays as a rule must be 
one of the latest operations of the whole device production cycle because otherwise the QD arrays would 
be destroyed by further high-temperature anneahngs. High-temperature processes during Ge/Si 
heterostructure formation on the late phase of chip production would, in turn, certainly wreck a circuit 
already formed on the crystal. Therefore, lowering of the array formation temperature down to the values 
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of < 450°C, as well as decreasing of the wafer annealing temperatures and times during the clean Si (001) 
surface preparation, is strongly required [20]. We refer to the Ge QD arrays and heterostructures based on 
them which satisfy this requirement as CMOS-compatible. 

Si (001) hydrogenation as a promising way of reduction of the surface cleaning temperature 

Development of a procedure of clean Si (001) surface preparation at lowered temperatures and/or by short 
thermal treatments is a keystone of creation of a CMOS compatible process of nanoelectronic VLSI 
fabrication [18l[32]. One of the ways of solving this problem is surface hydrogenation during wet chemical 
etching with subsequent hydrogen desorption from the surface in UHV ambient [27l[53]. In this connection, 
an issue of surface structure after these treatments becomes a task of primary importance taking into 
account a possible effect of Si surface atomic-scale roughness on formation of nanostructured elements 
(e.g., self-assembled Ge quantum dot nucleation on wetting layer in Ge/Si(001) 

heterostructures [20l[32l[54l[55] ) . In this section, we present data of our recent investigations conducted by 
means of STM and RHEED on preparation of clean Si(OOl) surfaces by hydrogenation and thermal 
desorption of hydrogen in an UHV MBE chamber after wet chemical etching by the RCA process [3T] . 

It is known [53l[56] that the temperature of surface cleaning depends on composition of etchants used for 
hydrogenation. Solutions based on HF, with pH varied from 2 to 7, are typically used for surface 
hydrogenation. A number of silicon hydrides form on the Si(OOl) surface by the reaction of hydrogen with 
Si, and the most typical ones are monohydride and dihydride (Figure [7j). A fraction of dihydride on the 
surface grows with the increase of pH; monohydride desorbs from the surface at higher temperature [53] . 

The main results of our studies are as follows. 

Explorations of hydrogenated surfaces [31 have evidenced that regardless of the type of solution used for 
surface hydrogenation, RHEED patterns correspond with unreconstructed 1x1 surface (Figure [8]). Broad 
streaks with pronounced 3D-related structure form the RHEED patterns for the samples etched in HF 
solutions (Figure [8^,b); high intensity of the Kikuchi lines indicates that the surface is highly smooth and 
ordered on macroscopic scale. Visible local enhancement of signal of the RHEED patterns takes place 
owing to overlapping of Kikuchi lines. The shapes of the streaks corresponding to the surface well 
developed on the monoatomic scale (3D spots) are detected in the patterns of the samples treated in NH4F 
solutions (Figure [8]i,e). According to STM data, the surface is more rough in the case of etching in the 
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NH4F solutions than in the case of hydrogenation in solutions based on HF. 

STM measurements have shown that clean Si (001) surfaces may be obtained as a result of hydrogen 
thermal desorption in the interval from 470 to 650°C ^J, but their roughness depends on chemical 
treatment applied for hydrogenation and temperature of subsequent annealing (Figures [9l and [TQ|) . If dilute 
HF solutions and annealings of moderate duration at temperatures higher than 600°C are applied, smooth 
surfaces with monoatomic steps and wide terraces are obtained (Figure [9^- d). The c(4 x 4) reconstruction 
is observed for such samples. If the duration of annealing at the temperature higher than 600°C is 
increased, SiC islanding may occur on the surface. Lower temperature of annealing gives rise to formation 
of a rough (2 x l)-reconstructed surface (Figure [9^-j). 

Application of solutions based on NH4F followed by any low-temperature annealing enables obtaining of 
clean rough Si(OOl) surfaces composed by narrow and short monoatomic steps (Figure [TQ|). The 
(2 X l)-reconstructed surface form as a result of annealing at the temperatures higher than 600°C 
(Figure fTOtl.e). 1x1 surface was observed after treatments at lower temperatures (Figure [10^, h). We 
would like to emphasize that annealings at the temperatures from 470 to 600°C result in formation of rough 
surfaces regardless of the applied chemical treatment (compare Figures [9t-j andfTOT-h): application of 
solutions containing NH4F always results in formation of more rough surfaces in comparison with surfaces 
of specimens treated in dilute HF solutions (Figure [TQ|). 

Notice that the (2x1) RHEED patterns were observed for the hydrogenated surfaces after annealing at 
800°C for 5 minutes and quenching |28] which were used as the reference samples with known surface 
structure. 

It should be noted also that comparison of the above STM and RHEED data makes one infer that RHEED 
cannot be applied as the only method of monitoring of the surface cleaning grade and the state of 
dehydrated surfaces [31j. RHEED patterns on the hydrated surfaces corresponded to the 1x1 structure. 
Surfaces cleaned as a result of subsequent annealings in the temperature interval from 470 to 650°C were 
(1 X 1) or either (2 x 1) or c(4 x 4)-reconstructed. Hence, in some cases, a type of the RHEED pattern did 
not change after thermal desorption of hydrogen, however forms of the patterns, which corresponded to the 
1x1 structure, were different before and after annealings. 

As of now, we suppose that clean Si(OOl) surfaces applicable for MBE formation of Ge/Si(001) 
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heterostructures can be obtained at low temperatures (as low as 470°C). However, it is not excluded that 
further lowering of temperature of the clean Si(OOl) surface preparation is possible, perhaps down to the 
temperatures as low as 400°C [53j. In the latter case. Si weak flux and formation of a buffer layer may be 
useful to prepare a perfect enough Si(OOl) surface before Ge deposition. 

Concluding this section, let us briefly consider the morphological peculiarities of the c(4 x 4)-reconstructed 
surface shown in Figure [9^, b. Figure [TTI presents magnified STM images of this surface. A structure 
observed in the images represents a mixture of a-c(4 x 4) and /3-c(4 x 4) modifications [57H58j: (2x1) and 
c(4 X 4)-reconstructed domains coexist on the surface (Figure [TTb): location of dimers forming the c(4 x 4) 
structure with respect to the dimers of the (2x1) structure is also seen; ad-dimers in both epitaxial and 
non-epitaxial orientations are seen in (Figure [TTb). The /3-c(4 x 4) modification prevails on the surface 
shown in Figure [TTti which is only partially occupied by c(4 x 4). It is seen that the presented data are in 
good agreement with the model of the c(4 x 4) structure proposed by Uhrberg etal. |57[ [58]. 

Conclusions 

Is 600°C a fundamental value of temperature for Si and Ge (001) surfaces? 

Concluding the article we would like to draw the reader's attention to the fact that many of the processes 
described above or in the cited articles have some critical temperature close to 600°C. Thus, the phase 
transition between 2x1 and c(8 x 8) reconstructions occurs around this temperature [T9"^"28J. Exploration 
of dehydrogenation of the Si:H samples shows that clean surfaces obtained by annealing at the 
temperatures > 600°C are formed by wide terraces with monoatomic steps; the c(4 x 4) reconstruction 
appears at these temperatures [57J. Annealings at the temperatures < 600°C result in formation of rough 
surfaces composed by narrow and short steps. Ge QD arrays deposited by MBE at the temperatures < and 
> 600°C also strongly differ in both cluster compositon and nucleation. Bimodal hut arrays form at low 
temperatures, whereas arrays grown at high temperatures are composed by pyramids and domes. The 
low-temperature clusters nucleate by formation of strictly determined 2D structures composed by dimer 
pairs and longer chains [20^32j. There are two alternative scenarios of cluster formation at high 
temperatures: (i) similarly to the low-temperature nucleation of pyramids and (ii) by {105}-faceting of the 
Ge shapeless heaps. An assumption arises from these examples that the temperatures do not coincide 
accidently, but some changes happen in the processes of migration of Si and Ge adatoms over the (001) 
surface around 600°C. 
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Summary 

In summary, issues of morphology, nucleation and growth of Ge cluster arrays deposited by ultrahigh 
vacuum molecular beam epitaxy on the Si (001) surface are considered in the article. Difference in 
nucleation of quantum dots during Ge deposition at low (< 600°C) and high (> 600°C) temperatures is 
studied by high resolution scanning tunneling microscopy. The atomic models of growth of both species of 
Ge huts — pyramids and wedges — are proposed. The growth cycle of Ge QD arrays at low temperatures is 
explored. A problem of lowering of the array formation temperature is discussed with the focus on CMOS 
compatibility of the entire process; a special attention is paid upon approaches to reduction of treatment 
temperature during the Si (001) surface pre- growth cleaning, which is at once a key and the 
highest-temperature phase of the Ge/Si(001) quantum dot dense array formation process. The 
temperature of the Si clean surface preparation, the final high-temperature step of which is, as a rule, 
carried out directly in the MBE chamber just before the structure deposition, determines the compatibility 
of formation process of Ge-QD-array based devices with the CMOS manufacturing cycle. Silicon surface 
hydrogenation at the final stage of its wet chemical etching during the preliminary cleaning is proposed as 
a possible way of efficient reduction of the Si wafer pre-growth annealing temperature. 
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Endnotes 

^The ideas underwritten in this paragraph have already been proposed by us [59] and A. V. Dvurechenskii. 
^In addition, an analytical chamber equipped with secondary ion mass spectroscopy (SIMS), Auger 
electron spectroscopy (AES) and X-ray photoelectron spectroscopy (XPS) is also available in the 
instrument. A Knudsen effusion cells for layer doping by boron during deposition is installed in the MBE 
chamber but it was not used in the described experiments. 

*^We express our acknowledgement to the anonymous colleague who proposed this explanation. 
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Figures 




Figure 1: 

Figure [1]- STM empty-state images of hut nuclei on Ge WL formed at different temperatures: 

(a)-(c) pyramid (1) and wedge (2) nuclei on the adjacent M x N patches of WL; Tgr = 360°C, hce = 6 A; 
the structural models [20l[32] are superimposed on the corresponding images in (a); (b),(c) pyramid nuclei 
on WL formed at low temperature (Tgr = 360°C): (b) hce = 5.4 A; (c) hce = 6 A; (d),(e) a pyramid nucleus 
on WL formed at high temperature, hce = 5 A: (d) Tgr = 600°C, 43 x 37 nm; (e) Tgr = 650°C, 7.8 x 6nm. 



20 



Figure 2: 



Figure [2]- Formation of {105} facets from shapeless areas with excess of Ge: 

STM images of different phases of faceting, Tgr = 650°C, hce = 5 A; (a) a shapeless Ge 'heap' without 
faceting, 150 x 141 nm; (b) at the outset of faceting, 64 x 64 nm; (c)-(d) after growth stoppage and 
anneahng at the growth temperature; (c) 72 x 72 nm; (d) 46 x 46 nm; (e) 23 x 23 nm. 
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Figure 3: 



Figure [3]- STM images of Ge WL on Si(OOl) at the outset of QD array formation: 

Tgr = 360°C, (a) hce = 4.4 A, = —1.86 V, It = 100 pA, neither hut clusters nor their nuclei are 
observed; (b) hce = 5.1 A, Us = +1.73 V, h = 150 pA; (c) Us = +1.80 V, h = 100 pA; (d) Us = +2.00 V, 
It = 100 pA. Examples of characteristic features are numbered as follows: nuclei of pyramids (1) and 
wedges (2) [1 ML high over WL patchs, Figure [I], [20l[32] smah pyramids (3) and wedges (4) [2 ML high 
over WL patchs] [I8l[20l[32lll9] . 
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Figure 4: 



Figured]- STM images of Ge WL on Si(OOl): 

Tgr = 360°C, the ordinary c(4 x 2) (c) and p{2 x 2) (p) reconstructions within the M x N patches are often 
observed simultaneously, (a) hce = 4.4 A, Us = —1.86 V, It = 100 pA, only the c(4 x 2) structure is 
resolved; (b) Hqq = 5.1 A, [/g = —3.78 V, It = 100 pA, both c(4 x 2) and p(2 x 2) structures are revealed as 
well as nuclei of a pyramid (1) and a wedge (2); (c),(d) hce = 6.0 A, Us = +1.80 V, It = 80 pA, both 
c(4 X 2) and p{2 x 2) reconstructions are well resolved; (e) hce = 5.1 A, [/g = —3.78 V, It = 100 pA, a 
pyramid nucleus on the c(4 x 2) reconstructed patch with the adjacent p{2 x 2) reconstructed patch. 
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Figure 5: 



Figure [5]- Models of Ge hut growth: 

(I) uniform addition of Ge atoms to four facets; (II) nun-uniform addition of Ge atoms to facets. 
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Figure 6: 



Figure [6]- STM images of a Ge QD dense array at different phases of its evolution from patched WL 
to 2D nanocrystalline layer: 

Tgr = 360°C, hce = (a) 4.4 A, before nucleation (see also Figure [3] for details of array nucleation at 

^Ge = 5.1 A); (b) 6 A, growing small huts, nucleation goes on; (c) 8 A, maximum density 

(- 6 X lO^^cm-^); (d) 10 A, maximum uniformity, large huts start to coalesce; (e),(f) 14 A, huts go on 

coalescing; (g) 15 A, 2D layer starts to form; (h),(i),(j),(k) 18 A, 2D nanocrystallyne film grows, chaos of 

faceted hillocks and pits (i) is observed; however, Ge WL (M x A/')-patched structure is clearly resolved on 

bottom of pits (j),(k). 
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Figure 7: 

Figure [7]- Structure of a hydrated Si(OOl) surface: 

(a) mono- hydride and (b) di- hydride. 
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Figure 8: 



Figure [8]- RHEED patterns and STM images of Si:H surfaces obtained as a result of different 
chemical treatments: 

(a)-(c) after hydrogenation in dilute HF; (d)-(f) after hydrogenation in buffered HF + NH4F; RHEED 
patterns: E = lOkeV, (a),(d) [110] azimuth, (b),(e) [010] azimuth; STM empty-state images: (c) 
100 X 100 nm, = +1.9 V, h = 100 pA; (f) 88 x 88 nm, = +2.0 V, h = 100 pA. 
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Figure 9: 



Figure [9]- STM images and RHEED patterns of Si:H surfaces obtained as a result of hydrogenation in 
dilute HF after different heat treatments: 

(a),(b),(e)-(h) STM empty-state images; (a) 650°C for 8min, 57 x 57nm; (b) 610°C for lOmin, 41 x 41 nm; 
(c),(d) corresponding RHEED patterns, E = lOkeV: (c) [110], (d) [010]; (e) 570°C for 20min, 
101 X 101 nm; (f) 550°C for 30min, 66 x 66 nm; (g) 530°C for 35min, 41 x 41 nm; (h) 500°C for 35min, 
49 X 49nm; (i),(j) corresponding RHEED patterns, E = lOkeV: (i) [110], (j) [010]. 
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Figure 10: 



Figure [To]- STM images and RHEED patterns of Si:H surfaces obtained as a result of hydrogenation 
in NH4F or HF + NH4F solution after different heat treatments: 

(a)-(c),(f) STM empty-state images; (a) NH4F, 650°C for 5min, 40 x 40 nm; (b) HF + NH4F, 610°C for 
lOmin, 56 x 56 nm; (c) NH4F, 610°C for lOmin, 88 x 87 nm; (d),(e) corresponding RHEED patterns, 
E = lOkeV: (d) [110], (e) [010]; (f) NH4F, 550°C for 35min, 60 x 60 nm; (g),(h) corresponding RHEED 
patterns, E = lOkeV: (g) [110], (h) [010]. 
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Figure 11: 



Figure [IT]- STM images of the Si(001)-c(4 x 4) surface: 

(a) empty states, 19 x 14 nm, = +2.5 V, It = 120 pA; (b) empty states, 19 x 12 nm, = +2.0 V, 
It = 120 pA; (c) filled states, 22 x 15 nm, = -3.9 V, It = 150 pA; (d) filled states, 19 x 20 nm. 
Us = —3.9 V, It = 150 pA. As it is clearly observed in (a), the structure is composed by a mixture of the 
a-c(4 X 4) and /3-c(4 x 4) modifications; it is seen in (c) that the c(4 x 4) and (2 x 1) reconstructions 
coexist on the surface; location of dimers forming the c(4 x 4) structure with respect to the dimers of the 
(2x1) structure is also seen; ad-dimers in both epitaxial and non-epitaxial orientations are seen in (c). 
The /3-c(4 X 4) modification prevails in (d) which is only partially occupied by c(4 x 4). 
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